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We have studied the 3-flavor, finite temperature, QCD phase transition with staggred fermions on an A''t = 4 
lattice. By studying a variety of quark masses we have located the critical point, rric, where the first order 3-flavor 
transition ends as lying in the region 0.32 < rric < 0.35 in lattice units. 



1. INTRODUCTION 

The flavor and mass dependence of finite tem- 
perature QCD phase transition has been studied 
extensively in the past [Q. Three flavor QCD is 
especially interesting because of the existence of 
a second order end point of the first order phase 
transition line which is conjectured to belong to 
the 3-d Ising universality class &. It is similar to 
systems such as the liquid-gas [|, the 3-d Ising 



((AM)^) 
((AM)-; 



^, in which M is the magnetic- like 



and the SU(2) 



|5|, the 3-d 3-state Potts model 
Higgs modelQ. It's worth mentioning that the 
field mixing phenomenon Q is quite important in 
explaining quantitatively some critical behaviors 
of these systems, but it does not have an impor- 
tant effect on the behaviors we are discussing. 

In this study, we extend our previous effort 
m to locate the end point using the standard 
Wilson gauge action and unimproved staggered 
fermion action. There are many approaches to 
locate the end point, which is charaterized by a 
diverging correlation length. In this paper, we 
study the discontinuity of the order parameter 
{'ipip) [H Etnd the meson screening masses (includ- 
ing the scalar singlet a meson) 1^ along the first 
order phase transition line and extrapolate to the 
end point. These two approaches works well in 
the region with relatively strong first order phase 
transition. Another method called the Binder in- 
tersection method ^ is based on the finite size 
scaling of symmetry sensitive quantities such as 
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observable. This method can be used to study 
the region close to the end point and does not re- 
quire a knowledge of critical exponents. However, 
it's rather expensive to obtain a large amount of 
statistically independent samples from full QCD 
simulations using existing evolution algorithms, 
especially in the critical region where the critical 
slowdown is problematic. The Binder intersec- 
tion method is currently under study and is not 
reported here. 



2. EVOLUTION 

To determine the first order phase transition 
line, we evolve the system starting from both cold 
and hot initial configurations. We use the hy- 
brid Monte Carlo R algorithm with a trajectory 
length of 0.5 molecular dynamic time and step 
sizes ranging from 0.01 to 0.0125 depending on 
quark mass. If a strong first order phase tran- 
sition exists, we should be able to observe two- 
state coexistence at large enough volume. This 
is indeed found at the small quark masses 0.015 
and 0.02 (Fig. 0) on an L = 16 lattice. At 
m = 0.025 , tunneling between the two phases 
happens in our simulation on an L = 16 lattice 
(Fig. H). However, this tunneling is not found 
on a larger, L = 32 lattice (Fig. ||), indicating 
that the phase transition is still first order. For 
m — 0.035, no two-state signal is found on either 
L = 16 (Fig. 11) or L = 32 lattices. The system 
also shows significant critical slowdown behavior 
and large fiuctuations, which indicate a long cor- 



relation length and closeness to a second order 
critical point. Crossover behavior is observed at 
TO = 0.04. 
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Figure 1. The (ipip) evolution and histogram for 
TO = 0.02 and L = 16 at /3c 5.1235 
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Figure 2. The {^pip) evolution and histogram for 
TO = 0.025 and i = 16 at /3c = 5.132 

3. SCREENING MASS 

Theoretical study ||] shows that the scalar sin- 
glet meson becomes massless while other mesons 
stay massive at the end point. To improve statis- 
tics, we measure the screening correlators in all 
three spacial directions and average them. The 
disconnected diagram needed for the flavor singlet 
case is measured using a noise estimation tech- 
nique [|l0|. Wc did a careful study of the effect 
of the type and number of noise sources to ensure 
that the noise from noise estimator is smaller than 
the gauge configuration noise. Figure || shows 
the time dependence of a typical correlator we 
obtained. The meson screening masses are mea- 
sured at the critical (3 value where a two-state 
signal is observed. It becomes difficult for large 
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Figure 3. The {ipi') evolution and histogram for 
TO = 0.025 and L = 32 at /3c = 5.132 
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Figure 4. The (iptp) evolution and histogram for 
TO = 0.035 and L = 16 at /3c = 5.151 

quark masses, near TOc which needs very large 
volume to reduce frequency of tunneling events. 
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Figure 5. Screening correlators 

Figure |6| shows the result of the pseudoscalar 
and scalar flavor singlet (a) screening masses. 

Extrapolating using the form Ma- oc (toc — 
to)^/^, we get Toc = 0.034(3) from the conflned 
phase and TOc — 0.032(2) from the deconfined 
phase. 
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Figure 6. (Mscrccning)^ VS. quark mass. 

4. i^i;) DISCONTINUITY 

The discontinuity in the order parameter {ijjip) 
decreases with increasing quark mass and finally 
disappears at rric. We fit the data with the mean 
field exponent (3 = 1/2 (Fig. defined as 

A(V;V') oc (toc - m)'^, and find = 0.0334(17). 
Fitting to the 3-d Ising exponent of 0.3285 gives 
a much poorer x^- 
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Figure 7. Discontinuity in (ipip) with a mean 
field extrapolation. 



5. SCALE SETTING 



In conclusion, we have obtained the approxi- 
mate location of the end point by studying the be- 
havior of order parameter discontinuity and me- 
son screening masses as quark mass approaches 
the end point of first order phase transition line. 
From our scale setting calculation, we conclude 
that there is no first order phase transition for 
three flavor QCD at the strange quark mass. 
Higher statistics and larger volume studies in the 
critical region are needed to determine more ac- 
curately the critical point, critical exponents and 
universality class. 
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We have done a scale setting calculation on an 
8'^ X 32 lattice. The pseudoscalar to vector mass 
ratio is 0.3134(34) and 0.3666(32) for m = 0.025 
and m = 0.035, respectively. If we use the phys- 
ical p mass for the vector meson, this implies an 
SU(3) symmetric pseudoscalar mass of 281 Mev 
at m = 0.035. 



